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EFFECTS OF VEGETATION CANOPY 


STRUCTURE ON REMOTELY SENSED CANOPY TEMPERATURES 

D. S. Kimes 
ABSTRACT 

Remote sensing of vegetation temperatures is a promising technique for inferring plant water 
stress and yield on a large spacial scale. The effects of vegetation canopy structure on thermal 
infrared sensor response need to be understood before vegetation surface temperatures of canopies 
witli low percent ground cover can be accurately inferred. The response of a sensor is a function 
of vegetation geometric structure, the vertical surface temperature distribution of the canopy 
components, and sensor view angle. Large deviations between the nadir sensor effective radiant 
temperature (ERT) and vegetation ERT for a soybean canopy were observed ^.hroughout the 
growing season. Tlie nadir sensor ERT of a soybean canopy witii 35% ground cover deviated from 
the vegetation ERT by as mucii as 1 1°C during the mid-day. These deviations were quantitatively 
explained as a function of canopy structure and soil temperature. Remote sensing techniques 
wiiicli uniquely determine the vegetation canopy temperature(s) from the sensor response need 
to be studied. 
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EFFECTS OF VEGETATION CANOPY 
STRUCl'URE ON REMOTELY SENSED CANOPY TEMPERATURES 

INTRODUCTION 

Advances in thermal scanner technology in the last few years Iiave created great interest in 
utilizing the tJiermal infrared region to gain additional information concerning tlie status of earth 
resources. Two NASA satellites have been launched which contain thermal scanners: Landsat-C 
and the Heat Capacity Mapping Mission (HCMM). Future missions with thermal IR scanners are 
being considered. Currently there is great interest in utilizing the thermal infrared region to make 
inferences about vegetation canopy characteristics. 

There have been a number of studies concerning the use of vegetation surface temperatures 
along with other variables to infer the water status of the vegetation canopies (Blad and Rosenburg, 
1976;PIeilman et al., 1976; Jackson et ah, 1977; and Ehrlcr et ah, 1978). Jackson etah, 1977 cited 
other related studies. In addition there has been some work using soil surface temperature to predict 
the water status of bine soil as reviewed by Idso et al. (1975). Both temperature and moisture of 
t];e canopy components (leaves and soil) are of primary importance in determinuig crop yields. As 
a specific example Idso et ah (1977) used wheat canopy temperatures and auxiliary air temperature 
measurements as an index to plant water stress for the period from head emergence to the cessation 
of head growth. This index was successfully used to estimate crop yield. The study demonstrated 
a promising empirical relationship to schedule irrigations and predict yield by remote sensing 
techniques. 

There is a need for monitoring the temperature of vegetation foliage and soil temperatures 
during the early stages of growth before tlie vegetation completely covers the ground. Such infor- 
mation is useful in predicting the maximum potential yield of a crop. For example, Idso et ah 
(1977) hypothesized that information on wheat canopy temperatures and soil surface temperatures 
during early crop stages would be useful in predicting the maximum green leaf area izidex which 
may detennine the final potential grain yield. Canopy temperatures early in the growing season 
would also be useful in scheduling tlie first few irrigations in arid lands. 
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The effects of the canopy’s geometric stnicturc on tlie sensor response need to be understood 
and corrected for before empirical relationships dependent on vegetation temperatures such as 
presented above can be accurately applied to vegetation canopies with relatively low percent ground 
cover. In this study the theoretical aspects of these effects were explored and tlie magnitude of 
the variability of some of these effects were documented for a soybean canopy, 

BACKGROUND 

The geometric structure of a vegetation canopy can be mathemetically described by such 
physical characteristics as the distribution and density of plants on the ground, foliage-area-index, 
and foliage angle frequency distribution as discussed by deWit (1965). From these parameters the 
probability of gap (PGAP) through horizontally infinite canopy layers as a function of view angle 
can be estimated (Nilson, 1971). 

The PGAP function of a vegetation canopy is an important parameter in determining tliennal 
infrared, radiant transfers. In the thermal infrared region the emissivity and absorptance values of 
natural objects approach 1,0 (Idso et ah, 1969). Thus for first order approximations, radiant trans- 
fers within a vegetation canopy (including the ground) are simplified in that reflections are minimal 
(Ross, 1976). As a consequence, only direct line emissions from the source to the sensor (e.g. 
a leaf surface within the canopy having an unobstnicted path to the sensor) need to be considered. 
The vertical surface area distribution of canopy components (leaves and soil) which has a direct 
line emission path to the sensor can be conveniently described by the PGAP function as presented 
below. 

The contribution of thermal radiance from tlie soil and vegetation canopy components to the 
sensor ii- a function of canopy geometry, spatial distribution of canopy component surface tem- 
pertures, and sensor view angle. Kimes et ah (J 979a) presented a numberical equation for esti- 
mating the thermal radiance of a vegetation canopy as a function of view angle. The canopy was 
abstracted into a number (n) of horizontally infinite layers (Figure 1). A generalized form of the 
equation is; 
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L(0,0) = 7r“‘ * S fen PGAPk (0,0)) • PHIT, (0,0) *e, • a -Tfl 
i=l Lk=0 J 

where 

L(0,0) = Thermal infrared radiance of a canopy as a function of view angle (0, 0) where 0 

and 0 are the inclinatior; and azimuth angles respectively (W*M“^ • SR"^ ) 
n - Number of discrete layers in the canopy system including the ground 

PGAP^(0,0) = Probability of gap in the direction 0, 0 for layer k. By convention 
PGAPo (0,0) = 1.0 for all 0 and 6. 

PHITj(0,0) = Probability of hit in tije direction 0,0 for layer i and PGAPj; (0,0) = 1-PHITj (0,0) 
for .i=k. By convention PHITn (0,0) = 1.0 for all 0 and 0 where 
n denotes the ground. 

6j = Mean component emisivity in layer i 

Tj = Mean component surface temperature in layer i (°K) 

c = Stefan-Boltzmann constant(W*M"^*°K‘’). 

The underlying assumptions and justification of the above formula are presented by Kimes et al. 
(1979a). 


^ M 

The (II PGAPjj (0,0))*PHITj(0,0) term can be thought as the proportion of projected leaf 
k=o 

area which is in the direct line of sight to the sensor for a particular view angle described by 0, 0 
and for a specific layer i. The remaining portion of the equation is just the Stefan-Boltzman 
equation for emitted exitance. Thus, Equation 1 sums the contribution of thermal infrared 
radiance from each layer of the canopy in a particular direction (0,0). The simplest case of Equation 
1 was used in the following analysis. 

The observed vertical canopy temperature distributions are a result of a number of simul- 
taneous energy transfers. These transfers include foliage transpiration, soil and foliage evaporation, 
soil and foliage solar absorption, thermal infrared emission and absorption for soil and foliage. 
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soil conduction, and soil and foliage convection. These energy transfer functions between tlie cn- 
viromnent and the vegetation canopy are discussed by Rosenberg (1974), Gates (1975), and Ross 
(1976). 

Many of these energy transfers are controlled by the vegetation canopy geonletr 3 ^ For 
example, the proportion of spectral solar irradiance which is absorbed by each canopy layer is 
controlled as a function of canopy geometry, optical properties of canopy components, and solar 
zenith angle (Kimes et al., 1979a). The variation of absorbed solar irradiance in canopy layers due 
to changing solar irradiaiice conditions is paraded by a change in layer temperature. These relation- 
ships, however, are complicated by other energy transfers (convection, transpiration, evaporation, 
and conduction) which occur simiilataneously as discussed by Gates (1975). These transfers are a 
complex function of wind speed, relative liumidity, air temperature, solar irradiance conditions, 
vegetation geometry, several water relation factors of the vegetation and soil, and the optical and 
thermal parameters of the vegetation and soil. 

The remote sensing problem is to relate the response of a thermal infrared sensor to the plant 
and/or soil conditions of interest. However, as discussed above there are a large number of environ- 
mental and botanical permutations wliich effect the sensor response and mask the desired target 
characteristics of interest. Presently, there are no studies which attempt to explore tlie compre- 
hensive relationsliips between the observed variability in sensor response to the underlying enviro* 
mental and botanical conditions. 

In this study the effects of canopy geometry on the response of a nadir looking sensor were 
documented throughout the growing season for a soybean canopy. The observed variability in the 
response of a nadir sensor was quantitatively described by applying the simplest case of Equation I . 

STUDY APPROACH 

Radiometric thermal measurements of a soybean canopy at the USDA Beltsville Agricultural 
Research Center, Beltsville, Maryland were taken during the 1978 field season. The data were 
collected fmd reported by UHL et al. (1979). Tlie autiiors presented supporting microclimatic and 
agronomic data. Effective radiant temperatures (ERT) were taken with a Barnes PRT-5 Precision 
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Radiation Thermometer (field of view = 2°, accuracy = 0 5’c). All measurements were taken at 
approximately 1330 hours (Eastern Daylight Time) foi 10 days throughout the growing season. 

For each measurement period canopy ERT, ground ERT, and composite sensor ERT were 
obtained. Canopy ERT was measured at four random points with the field of view (FOV) of the 
radiometer entirely on vegetation components. A mean canopy ERT (T^,) was calculated. Ground 
ERT was measured at four random points with only the soil surface in the FOV of the radiometer 
and the mean ground ERT (Tq) was calculated. A composite sensor ERT (Tg) of the canopy was 
measured from a 10 m high platfonn. The Tg measurement was taken in the nadir direction. The 
and Tq measurements were assumed to be the mean surface temperatures of the vegetation 
canopy and ground, respectively. The Tg measurement represented the composite response from 
the scene composed of vegetation and ground. 

The percent ground cover of the vegetation was estimated ocularly by tlie same individual 
throughout the growing season. These percent cover measurem-jnts were used as an estimate of 
probability of hit (PHIT) in the nadir direction. PGAP and PHIT values are related as PGAP = 
l-PHIT. 

The variability observed in the Tg measurernent were explained quantitatively as follows. 
The canopy was abstracted as a simple two layered system (vegetation and soil). Only a nadir 
view angle of the sensor was considered. Equation 1 was reduced to tliis specific case and was 
expressed exclusively in terms of ERT ratlier than true surface temperatures and emissivities. The 
simplified equation is: 

T| = PGAP • T^ t (l-PGAP)Tj ^2) 

where 

A 

Tg = theoretical effective radiant temperature of a nadir looking sensor 
Tq = mean effective radiant temperature of ground 

T(, = mean effective radiant temperature of the canopy foliage 

PGAP = probability of gap through the canopy in the nadir direction. 
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The above measurements of Tq, T^,, and PGAP for tlie 10 measurement periods were used to 

A A 

calculate the respective values. Comparisons of the theoretical Tg and measured Tg were made, 
RESULTS AND DISCUSSION 

In several research applications as nored above, one would like to infer a mean canopy tem- 
perature from the Tg measurement. Figure 2 sli nvs the measured sensor ERT (Tg) versus the 
measured mean canopy ERT (T^,) as a function of percent cover. The dashed line in Figure 2 
and all subsequent figures is a reference line wiiere the dependent variable equals the independent 
variable. Equation 2 was used to explain the observed deviations between Tg and T^,. Equation 
2 was transformed as follows. 

T^ = rOAP • (T^ - tJ) + Tj . (3) 

From Equation 3 it is evident that only at relauvely hudi PGAP values and high (T^ - T^) differ- 
entials is the Tg significantly different than T^,. /'iS seen in Figure 2, the measured sensor ERT (Tg) 
significantly deviated from the canopy ERT (T^,) only for those measurements with high PGAP 
values (Low percent cover). The magnitude of tiiese deviations were directly related to the 
(Tq ~Tj^) differential as presented in Figure 3. 

The (Tq - Tj,) differential can vary widely. Canopy ERT (T^,) closely approximates air 
temperature during low water stress as reported by Gates (1966), Monleith (1973), Uhl et 
(1979), and others. The ground ERT (Tq), however, can deviate greatly frem air temperature. 
Uhl et al. (1979) reported that tlie surface temperature of bare soil surfaces was occasionally as 
much as 20°C higher than ambient air temperature. Tlie difference in the soil temperatures 
(Figure 3) can be qualitatively explained by variations in solar irradiance, soil moisture, canopy 
geometry, and wind speed. 

Solar irradiance absorbed by the grouted is strongly affected by geometric and optical pro- 
perties of the overlying vegetation canopy. The physical principles were discussed by Oliver and 
Smith (1974) and Ross (1976). From simulated results and literature studies, Kimes et al. (1979b) 
showed for photosynthetically active wavelengths that spectral absorption of the ground under 
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vegetation canopies was higlily variable as a function of canopy »,'eoinetry, and solar zenith angle. 
This variability in ground solar absorption changes the soil energy budget which effects the soil 
surface temperature. Uhl et al, (1979) reported that summer sell temperatures were generally 
reduced several degrees centigrade by the presence of vegetation. At the same time, however, the 
authors report a decrease in wind speed near the ground within vegetation canopies as opposed to 
bare ground. Thus, convectional transfers for the canopy and the ground varies as a function of the 
canopy geometry and environmental conditions, Uhl et al. (1979) also reported that soil moisture 
greatly influenced the soil surface temperature. The energy dynamics of the soil were discussed 
further by Geiger (1965) and Rosenberg (1974). 

A ^ 

Figure 4 represents the measured sensor ERT (Tg) versus the theoretical sensor ERT (Tg). 
Most of the variability in Tg was explained by Equation 2. The root mean square of the deviations 

A 

between Tg andTg was 1.4 C. 

Thcjiiean canopy temperature (T^,) is the variable \Vhich is often desirable to infer from tiie 
composite sensor ERT (Tg), As demonstrated above soil temperature can have significant effects 
on the composite sensor signal. Remote sensing can readily obtain a nadir looking Tg . However, 
to uniquely detennineT^, other information is required. Research is needed to e-xplore possible 
approaches to remotely correct for this composite effect. One possible approach Is to use additional 
infonnation gained by off nadir view angles. Such a technique is suggested by Equation 1. These 
possibilities vill be explored in a future article. 

The above analysis and discussion was based on the assumption that one mean canopy tem- 
perature (T^) is characteristic of the entire canopy. However, vertical temperature gradients do 
occur in vegetation canopies. Bauer et al. (1977) presented the vertical temperature profiles of two 
wheat canopies for various measurement periods. During the afternoon mean foliage temperature 
differentials (top-bottom layer) were as much as ~1 1°C. The authors suggested that such infor- 
mation could be used to improve the accuracy of water relation models. There is the possibility 
of remotely sensing mean canopy layer temperatures using a series of off nadir view angles 
(Equation I), 
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Not only do canopy temperature gradients occur but also air temperature gradients within 
the canopy. Bergen (1971) showed that air temperature differentials within lodgepole pine canopies 
can be as much as 4-5°C in the verskal profile for clear, sunny days. Bauer ct al. (1977) found 
tiiat air temperatur<. differentials within wheat canopies can be as higli as 7.5°C. These observed 
variations have important implications in developing and applying relationships dependent on air 
temperature. 

SUMMARY AND CONCLUSIONS 

Past studies have utilized vegetation surface temperatures to infer the water status and yield 
of vegetation canopies. Remote sensing of canopy temperatures is a promising technique for moni- 
toring plant characteristics on a large spatial scale. The tliermal infrared sensor response of a vege- 
tation canopy is not a simple measurement cT foliage temperature but is a function of canopy 
geometric structure, vertical foliage surface temperature distribution, soil temperature, and sensor 
view angle. The smface temperature distribution of foliage is a complex function of a number of 
simultaneous energy transfers. Many of these transfers are dependent on canopy geometric structure. 

The soybean data and analysis documented the effects of the geometric structure of a 
simple two layered system (ground and vegetation) on the response of a nadir sensor (Tg) througli- 
out the growing season. For vegetation canopies with relatively low percent cover, the ground 
surface temperature (T^ ) can significantly influence the sensor response (Tg). For a soybean canopy 
with 35% ground cover Tg exceeded the canopy temperature (T^,) by 1 1°C. T^. can be highly 
variable due to variations in absorbed solar irradiiince and soil moisture. Soil temperatures (T^) 
as high as 20®C above air temperature have been observed for bare soil. Tq for a soybean canopy 
with a 35% ground cover was as higlr as 15°C above canopy temperature. This solar heating 
plienomena is controlled by the vegetation geometry, solar zenith angle, optical properties of the 
soil surface and vegetation, and physical properties of the soil matrix. The deviations of Tg from 
T^ were quantitatively explained by variations in soybean canopy geometry and soil temperature. 

The root mean squared deviation between the calculated theoretical Tg and tlie measured Tg was 1 .4. 
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Remote sensine techniques which comect for or circumvent these geometric effects need 

to be deveioped before remote sensing of foiiage terapemtums can be eccurately infetred throughout 

the growing season. First order corrections are needed which separate the composite sensor signai 

into mean soii and vegetation temperatures. For mom deveioped canopies, second order eoirections 

which separate the sensor response into a vertical pronie of mean canopy layer temperatures would 
be desirable. 
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Figure I. Abstraction of a vegetation canopy into four layers. Any number 
of canopy layers containing canopy components (leaves, stems, 
reproductive structures, and soil) can be abstracted. 
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Fijiuro 2 . Measured sensor FRT (Tj.) versus the measured mean canopy FRT (T^) 
of a soybean canopy for the 10 measurement periods ( 1330 hours) 
throughout the gnawing season. The number for each data point 
is the corresponding percent cover. 
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Figure 3. Measured mean ground ERT (T^ ) versus the measured mean 
canopy ERT (Tj.) of a soybean canopy for the 10 measurement 
periods ( 1330 hours) throughout the growing season. The number 
for each data point is the corresponding percent cover. 
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Figim* 4. Moasurcd sensor I R 1 ( I\,) versus the theoretieal sensr)r l-RT (Tj.) 
for a st^yhean canopy for the 10 measurement periods! 1330 liours) 
throughout the growing season. The number for each data point 
IS the eorrespoiuling percent cover. 
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